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On the basis of an electronic model with separable attractive interaction, the precursors at high 
temperature and strong coupling of the d-wave superconducting state are investigated in the one- 
particle spectral function A(k,uj) and the total density of states p{u)), with the use of the self- 
consistent t-matrix approximation. In the low-density region, it is found that a gap-like structure at 
the Fermi level appears in A(k,Lj) and p{uj) above the superconducting transition temperature. It 
is shown that the pseudogap energy scale is determined by the binding energy of the Cooper-pair. 
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I. INTRODUCTION 

In the underdoped high-T c superconductors (HTSC), 
pseudogap (PG) behavior has rbeen widely observed 
in experiments, such as NMR,EI specific heatj3 and 
photoemission.fj All these phenomena can be basically 
understood as caused by the suppression of low-energy 
spectral weight in the temperature range T c < T < Apq, 
where T c is the superconducting transition temperature 
and Apg is a characteristic energy scale for the PG 
formation. This occurs both in the spin- and charge- 
excitation spectra. As a consequence, the problem is re- 
duced to the clarification of the origin of this spectral- 
weight suppression, namely, the physical origin of Apg- 

The angle-resolved photoemission (ARPES) spectrum, 
which is sensitive to the momentum dependence of the 
PG, has revealed that the PG phenomenon itself exhibits 
a ci-wave symmetry which iS|-smoothly connected to the 
d-w&ve superconducting gap.0 Moreover, the locus of the 
minimum gap position in momentum space traces the 
shape of the Fermi surface. From these results, it can be 
inferred that the energy scale for PG formation is closely 
related to the superconducting correlation. Then, one of 
the possible explanations for the PG behavior involves 
the discussion of possible "precursors" of the Cooper- 
pair formation above T c . Certainly there are other pos- 
sible scenarios .-that also lead to PG formation sugh as 
spinon-pairing,u antifexromagnetic spin fluctuation,!?! and 
fermion-boson model,EI but in this paper the focus will 
be precisely on the development of a PG at strong cou- 
pling due to the formation of electron bound-states at a 
temperature scale larger than the one corresponding to 
long-range superconducting pairing. 

Along this scenario, much effort has been_devoted to 
the investigation of the PG phenomena.0LJ However, 
there are few results in the literature leading to PG with 
d-wave symmetry, while the PG in the s-wave supercon- 
ductor has been intensively investigated on the basis of 
the negative-?/ Hubbard model. The popularity of the 
s-wave calculations as opposed to the more realistic d- 
wave case is mainly due to technical issues. The quantum 



Monte Carlo (QMC) simulation provides accurate infor- 
mation on the negative-?/ Hubbard model and with these 
results the validity of other diagrammatic method such 
as the self-consistent i-matrix approximation (SCTMA) 
can be checked. However, for the model with c£-wave at- 
traction (or the nearest-neighbor attraction), QMC cal- 
culations are difficult mainly due to sign-problems in the 
simulations and also because phase separation could oc- 
cur for a model with an attractive potential that acts at 
finite distances, contrary to what occurs in the attractive 
Hubbard model where the attraction is only on-siteB3 

In spite of these potential difficulties, here the d-wave 
PG is studied in order to contribute to the investigation 
of the energy scale Apg in HTSC. For this purpose, here 
an effective model with <i-wave separable attraction is an- 
alyzed, focussing our efforts into the low-density regime, 
for the following reasons. First, from a physical point of 
view, the low carrier density region is important because 
the underdoped HTSC regime as a first approximation 
can be described as a low-density gas of holes in an an- 
tiferromagnetic background. Previous numerical studies 
have shown that holes in such an environment behave 
like quasiparticles with a bandwidth renormalized to be 
of order J, the Heisenberg exchange couplingo Second, 
now from a technical viewpoint, it is known .that the 
SCTMA gives reliable results in the dilute limitO Then, 
the behavior of the spectral function can be safely inves- 
tigated in the low-density region. For these reasons in the 
present paper the average electron filling will actually be 
at most 10%. Note that our "electrons" below will sim- 
ply represent fermions interacting through an attractive 
potential, and thus they can be thought of as "holes" in 
the context of HTSC. 

As mentioned before, the preformed s-wave pairing fea- 
tures in the negative-?/ Hubbard model have been widely 
investigated as a prototype for PG formation in the un- 
derdoped HTSC. Besides the technical aspects already 
discussed, this seems to be based on the assumption that 
the difference s vs d in the pairing symmetry does not 
play an essential role in the PG formation. This may 
seem correct by observing the gap-like structure in the to- 
tal density of states (TDOS), because it appears around 
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the Fermi level irrespective of the pairing symmetry, al- 
though the actual detailed shape is different. However, 
recalling that the main features for the PG formation in 
the underdoped HTSC have been revealed using ARPES 
technique, the structure in the individual one-particle 
spectral function should play a crucial role. In fact, 
important differences between s- and <i-wave symmetry 
fairly clearly appear in the spectral function described 
below in our study. 

In this paper, it is reported that the Apg scale agrees 
with the binding energy of the Cooper-pair irrespective 
of the pairing symmetry. The main difference between 
s- and d-wave symmetry, appears in the momenta of 
preformed-pair electrons, K and — K. For the s-wave 
symmetry, K is always determined by the band structure. 
Namely, in the dilute limit, it is given by the momentum 
at the bottom of the band, k*. Since the attraction is 
uniform in momentum space, K is determined only by 
the kinetic energy for the s-wave case. On the other 
hand, for a strong attraction with d-wave symmetry, K 
is not given by k*, but is located at (ir, 0) and (0,7r), 
because the attraction becomes maximum at those mo- 
menta. Such a competition between the band structure 
and the strong attractive interaction leads to interesting 
features in the d-wave PG, while the s-wave PG simply 
follows the band structure. 

This paper is organized as follows. In section II, a gen- 
eral formalism to calculate the electronic self-energy in 
the SCTMA on the real-frequency axis is present. For 
the investigation of the PG structure for the d-wave at- 
traction, a technical trick called "the s-d conversion" is 
introduced. Section III is devoted to the results obtained 
with the formalism of Sec. II. Two types of band struc- 
tures are considered with k* = (0,0) and (7r,0), respec- 
tively. In section IV, the results are discussed. Finally 
in section V, after providing some comments, the main 
results of this paper are summarized. Throughout this 
paper, units such that h = ks = 1 are used. 



II. FORMULATION 

A. Hamiltonian 

Let us consider a simple model in which electrons are 
coupled with each other through a separable attractive 
interaction. The symmetry of the electron pair is con- 
tained in the attractive term of the model, but it is not 
necessary to write it explicitly in most of the formulation 
of this section, although it will become important for the 
discussion on the PG. The model Hamiltonian is written 
as 

k,k',q 



where Ck CT is the annihilation operator for an electron 
with momentum k and spin a, £k the one-electron energy, 
fj, the chemical potential, and Vk,k' the pair-interaction 
between electrons. The electron dispersion is expressed 
as 

£k = — 2t(cos k x + cos k y ) — At' cos k x cos k y , (2.2) 

where t and t' are the nearest and next-nearest neighbor 
hopping amplitudes, respectively. The pair-interaction is 
written as 

Vk,k' = -KfkJv, (2.3) 

where /k is the form factor characterizing the symme- 
try of the Cooper-pair. Note that a positive value of V 
denotes an attractive interaction throughout this paper. 



B. Self-consistent t-matrix approximation 

Now let us calculate the spectral function using the 
SCTMA. Since this method becomes exact in the two- 
particle problem, it is expected to give a reliable result 
in the low-density region.Ej In fact, this expectation has 
been already checked in the attractive Hubbard model by. 
comparing SCTMA results against QMC simulations E3 
Therefore the reliability of the SCTMA may also be ex- 
pected for the non-s-wave attractive interaction, even 
though the direct comparison with QMC results is quite 
difficult in this case. 
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FIG. 1. (a) Self-energy in the SCTMA. The hatched square 
and the solid line denote the f-matrix V and the renormalized 
Green's function G, respectively, (b) Diagrammatic represen- 
tation for F. The broken lines indicate the interaction. 

Consider first for completeness the imaginary-axis rep- 
resentation. In this formulation, the one-particle Green's 
function G is given by 

G(k, iu> n ) = . (2.4) 

iuj n - (e k - (i) - E(k, iu) n ) 

where w n — 7rT(2n + 1), n is an integer, and T the tem- 
perature. In the SCTMA, the self-energy E(k, iu) n ) is 
obtained with the use of the i-matrix given by the in- 
finite sum of particle-particle (p-p) ladder diagrams, as 
shown in Fig. 1. More explicitly, £ is expressed as 
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E(k,«j n ) = /^ T^^r(k + k',iw„ + w n O 

n' k' 

xG(k'>„,), (2.5) 
where T(q, zi^i) is the i-matrix, given by 

" 2 <j!>(q, iu m ) 



1 - V0(q,M/ m )' 



(2.6) 



Here i/ TO = 27rTm, with m an integer, and 0(q, W m ) is 
the p-p ladder, defined by 



(cbu/ m )=T££/i!G(k,iw n ) 

n k 

x G(q - k, iv m - iw n ). 



(2.7) 



Note that the Hartree term is neglected in the self-energy 
because it should be considered as included in the band 
structure. The Green's function can be calculated self- 



consistently using Eqs. (2.4)-(2.7). The chemical poten- 
tial is determined by 



l )/2 = T^^e^""G(k,^„) 



(2.8) 



where (n) is the average electron number density per site, 
and rj is an infinitesimal positive number. In order to 
obtain results on the real-frequency axis, Pade approxi- 
mants for the numerical analytic continuation from the 
imaginary-axis data are frequently used.EH However, in 
general, it is difficult to control the accuracy of the cal- 
culation by this procedure. 

In this paper, our efforts are focused on the direct cal- 
culation of the Green's function on the real-frequency 
axis.lla In this context, a self-consistent calculation for 
the spectral function 



A(k,w) = (-l/7r)ImG(k,w), 



(2.9) 



is carried out, where the retarded Green's function 
G(k, lo) is given by 



where p. v. means the principal-value integral. The t- 
matrix is 



r(q,w) = 



l-V^faw)' 
where Im(/>(q, uS) is given by 



(2.13) 



f to' 
Im <f>(q,uj) =^22 dw'fl tanh— A(k,u/) 

1, J 



x A(q_ — k, lj — a/), 



(2.14) 



and the real part of </>(q, to) is also obtained using the KK 
relation. The electron number is obtained through 

(n}/2 = ^ f dtoA(k,uj)f F (uj), (2.15) 



and the spectral function must satisfy the sum-rule 



1 = ^2 I du;A{k, w) 



(2.16) 



This will be a check for the accuracy of the numerical 
results presented here. 

In the actual calculation, the fast Fourier .transforma- 
tion is applied to accelerate the procedureEa The first 
Brillouin zone is divided into a 64 x 64 lattice and the 
frequency integration is replaced by a discrete sum in 
the range — 25t < to < 25t, dividing it into 512 small in- 
tervals. As a consequence, the energy resolution is about 
O.li, indicating the order of magnitude of the lowest tem- 
perature at which our calculations can be reliably carried 
out. When the relative difference between two successive 
iterations for A(k,uj) is less than 0.01 at each (k, to), the 
iteration loop is terminated. As for the sum-rule, it is sys- 
tematically found to be satisfied within 1%. This value 
is the limitation for the accuracy of the present calcula- 
tion, because the integral equation with a singular kernel 
is solved by replacing the integration procedure by a sim- 
ple discrete summation. 



G(k,w) 



w - (ek - m) - E(k,w) 



(2.10) 



The imaginary part of the retarded self-energy is ex- 
pressed as 

Im£(k,a;)=/^ f dw'[f F {uj') + / B (w + uj')} 

k' 

x A(k', u/)ImT(k + k',uj + uj'), (2.11) 

where f F (x) = \/{e x / T + 1) and f B (x) = l/(e x / T - 1). 
The real-part of E is obtained through the use of ImE in 
the Kramers-Kronig (KK) relation 

„ . f dto' Im E(k, uj') 
Re E(k,w = p.v. / K — — '-, (2.12) 

/ IT UJ — <jj' 



C. s-d conversion 

For a separable potential with d-wave symmetry, /k is 
given by 



/k = cos k x — cos k y 



(2.17) 



In this case, due to the prefactor fy, in Eq. ( |2.1l| ), E 
always vanishes along the lines k x — ±k y , leading to a 
delta-function contribution in the spectral function. In 
order to avoid this singularity, a self-consistent calcula- 
tion for the d-wave case was first attempted by imposing 
anti-periodic and periodic boundary conditions for the 
k x - and ^-directions, respectively. However, it was not 
always possible to obtain a converged self-consistent solu- 
tion in this case. Actually, it was quite difficult to control 
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such convergence even if the temperature T was slowly 
decreased from the high-temperature region in which a 
stable solution was obtained, or if the coupling V was 
adiabatically increased from the weak-coupling region. 

This difficulty is caused by the fact that E becomes 
negligibly small in the region around k x rs ±k y for the 
d-wave case, even if the strong-coupling value for V is 
set as high as V — 8t. If ImS becomes smaller that 
the energy resolution in the present calculation, which 
is about O.li, the sharp peak structure in the spectral 
function around k x « ±k y is not correctly included in 
the self-consistent calculation. This leads to a spurious 
violation of the sum-rule, indicating that technical prob- 
lems appear in reaching a physically meaningful solution 
for d-wave symmetry at low temperatures. 

In order to avoid this difficulty, a continuous change 
from s- to d-wave symmetry is here considered by intro- 
ducing a mixing parameter a such that 

fl = (1 - a) + a(cosk x - cosk y ) 2 . (2.18) 

Our calculations start at a = 0, i.e., for the pure s-wave 
case, in which a stable solution can be obtained easily 
in the SCTMA. Then, a is gradually increased such that 
the d-wave case is approached. If a physical quantity for 
the d-wave model is needed, an extrapolation is made by 
using the calculated results for the quantities of interest 
between < a < 1. 



III. RESULTS 

In this section, our results calculated with the use of 
the real- axis formalism are shown. Here the magnitude 
of the interaction V is fixed as V — 8t. 



a small complication that can be solved with more CPU 
and memory-intensive studies than reported here. 



0.02 



II 

- (a) 


1 1 ' 

' 9 
i ij 

i | 


I 

- 


V=8t 


i 




t'=0 


ij 




- T=0.5t 


ij 




<n>=0.02 


ij 
if 






i 






Jh 




~ 


















——i i i 1 


r— -i i 


i 



-2.0 -1.0 0.0 1.0 

to (units of t) 



0.10 



| 0.05 



0.00 





u 


1 1 1 1 1 1 


(b) 


ill 

Mil 


V=8t 




iii 


f=0 




ili 1 


T=0.5t 






<n>=0.02 




ill! 

Hi 1 


k=(0,0) 




1\ 
\\ 












1 -1 - - >' --J2S*a4. — i 



-4.0 -2.0 0.0 2.0 4.0 6.0 8.0 

to (units of t) 



A. Case of t' = 

Let us consider first the band structure with k* = 
(0,0). In Fig. 2(a), the total density of states p(u>) is 
shown, given by p(iv) — J2k ^4.(k, u>). The whole curve for 
TDOS is not presented in this figure, because its shape at 
a larger scale is quite similar to the non-interacting case. 
At a = 0, a gap-like feature at the Fermi level can be 
observed, although it is shallow. This result has been al- 
ready repotted in numerous previous papers using several 
tcchniques.c3 With the increase of a, the gap structure 
gradually becomes narrower and at the same time deeper. 
The TDOS extrapolated to a = 1 using the results for 
the as in the figure is not shown, because it becomes un- 
physically negative in some energy region. However, this 
is not a serious problem, because such a behavior is only 
an artifact due to the extrapolation using a small number 
of a-results and it will disappear if a approaches unity 
very slowly and calculations with higher-energy resolu- 
tion are performed. This problem is not present in the 
studies at t' = in the next subsection. Thus, this is 




-4.0 -2.0 0.0 2.0 4.0 6.0 8.0 
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FIG. 2. (a) Total density of states (TDOS) for several val- 
ues of a, the interpolating parameter between s— and d-wave 
models, in the case V = 8i, t' = 0, T = 0.5t, and (n) = 0.02. 
The solid, dotted, dashed, long-dashed, and dot-dashed 
curves denote the TDOS results for a = 0,0.2,0.4,0.6, and 
0.8, respectively, (b) The spectral function at k = (0,0). (c) 
The spectral function at k = (tt, 0). 
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In order to understand the observed changes in the PG 
behavior with the increase of a, special attention must 
be given to the spectral function A(k,oj). Let us first 
analyze the result at k = (0,0), shown in Fig. 2(b), in 
which two peaks are observed. The large peak above the 
Fermi level is due to the quasi-particle (QP) contribu- 
tion, because if the interaction is gradually decreased, it 
continuously changes into the expected non-interacting 
(5-function peak. Thus here it will be called "the QP 
peak" . However, note that another structure can be ob- 
served below the Fermi level, although it has only a small 
weight. As will be discussed in the next subsection, this 
originates from the peak structure in ImT. In this sense, 
it can be called "the resonant peak" due to the forma- 
tion of the bound pairJld When a is increased, the QP 
peak becomes sharper and the position of the resonant 
peak is shifted to the right side in Fig. 2(b), while the 
weight decreases. At a = 1, the resonant peak will likely 
disappear and only the 5- function QP peak will occur, 
since the self-energy vanishes due to the prefactor m 



it is here typically chosen as t' — —t but the results do 
not depend crucially on such a choice. 



Eq. ( p.llp . 

Although the weight for the resonant peak in A(k, ui) 
with k = (0, 0) decreases with the increase of a, it is actu- 
ally transfered to another A(h,u>) with k ^ (0,0). Then, 
let us next turn our attention to A(k, ui) with k = (tt, 0), 
shown in Fig. 2(c). In this QP peak is also ob- 

served, but the position is higher than that at k = (0, 0). 
The difference between the positions of those QP peaks 
is about At, namely equal to £( w ,o) — e (o,o)- It should be 
noted that another peak structure grows below the Fermi 
level with the increase of a. The position roughly agrees 
with the lower edge of the PG structure in the TDOS, 
suggesting that the PG structure for d-wave originates 
from A(k, uo) around k = (tt, 0). 

Let us summarize this subsection. Pseudogap features 
appear in the density of states both for s- and <i-wave 
models, but its origin is quite different. For the s-wave 
case, this structure is mainly due to the preformed pair 
of electrons around the point k = k* = (0, 0) at the bot- 
tom of the band. On the other hand, for the d-wave case, 
it originates from the pair of electrons at other k-points, 
especially, k = (tt, 0). In the case of strong attraction 
such as V — 8t, those electrons can exploit the effect 
of the attractive interaction, in spite of the loss of the 
kinetic energy. In other word, this difference is due to 
the competition between the kinetic and the interaction 
effects. 
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B. Case of t' ^ 

From the result for t' = 0, in order to obtain a large PG 
structure for d- wave symmetry, it is necessary to consider 
the band structure in which k*'s are located at (±7r,0) 
and (0, ±7r). The reason is that electrons around k = k* 
can exploit the kinetic as well as the pairing energy due 
to the strong attractive interaction. As for a value of t' , 



FIG. 3. (a) Total density of states (TDOS) for several 
values of a in the case of V = 8t, t' = — t, T — 2.0t, 
and (n) = 0.2. The solid, dotted, dashed, long-dashed, and 
dot-dashed curves denote the TDOS for a = 0, 0.2, 0.4, 0.6, 
and 0.8, respectively. The thick solid curve is an extrapo- 
lated to a = 1 result from the available TDOS's in the inter- 
val < a < 1. (b) The spectral function at k = (tt,0). (c) 
The spectral function at k = (0,0). 
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In the TDOS shown in Fig. 3(a), no structure around 
the Fermi level is observed for the s-wave case, but a peak 
appears below the Fermi level with the increase of a. It 
can be regarded as a sign of PG formation, but this inter- 
pretation becomes much clearer if AQt, to) is investigated. 
In Fig. 3(b), the change of A(k, lo) at k = k* = (tt, 0) is 
depicted when a is increased. In the pure s-wave case, a 
large QP peak can be observed, but it is difficult to find a 
resonant peak below the Fermi level. On the other hand, 
when the d-wave is approached by increasing a, the QP 
peak is gradually destroyed and the resonant peak grows 
strongly below the Fermi level. Then, the PG structure 
is much larger compared to that at t' = 0. Note that in 
this case, the extrapolation for the TDOS at a = 1 is 
quite successful, contrary to what occurs at t' — 0, be- 
cause the large size of the PG allows us to perform the 
calculation at a high temperature such as T — 2i, a sit- 
uation in which the structure in the TDOS is smoother 
than at t' = 0. 

For the case t' — —t, weight transfer in A(k, lo) is ob- 
served with the increase of a, but it occurs between the 
QP and the resonant peaks at k = (tt, 0). In order to 
confirm this idea, ^4(k, lo) at k = (0, 0) was studied as 
shown in Fig. 3(c). As expected, only the sharpening 
of the QP peak is observed as a is varied, because the 
strength of the attractive interaction at k = (0,0) be- 
comes weak with the increase of a. Note that a finite 
width for the QP remains at a. = 1, but it is only a nu- 
merical artifact. Actually at t' = —t, electrons around 
k = (0, 0) do not take part in the PG formation even 
for the s-wave case. Since electrons around k = k* gain 
both the kinetic and potential energies, the PG structure 
is determined only by those electrons. 



IV. DISCUSSION 
A. Energy scale for pseudogap 

From the results in the previous section, in addition 
to the QP peak, a resonant peak below the Fermi level 
in A(k*,w) has been observed, although its appearance 
depends on the value of t' and the symmetry of the pair 
interaction. These two peaks define the PG structure in 
A(]s*,lo) and also in the TDOS, although in the latter it 
is often difficult to observe due to the smearing effects of 
the sum over momentum of the individual one-particle 
spectral functions. Based on these observations, in this 
paper the PG energy Apq is defined by the width be- 
tween the QP and the resonant peaks in A(k*,Lo). Note 
that for t 1 — and a = I, the weight for the resonant 
peak in A(\<* , lo) will vanish, but in the limit of a — » I, its 
position approaches the lower peak of the PG structure. 

In order to elucidate the physical meaning of our 
Apo, the imaginary part of the self-energy is analyzed 
at k = k* , because its structure has a direct effect on the 
spectral function, given by 



A(k,w) 
I 



Im£(k, lo) 



- (e k - M) - ReS(k, lo)} 2 + [ImS(k, lo)} 2 



(4.1) 



For an intuitive explanation, it is not convenient to 
analyze the full self-consistent solution for ImS(k, lo). 
Rather the essential inf orma tion can be obtained by 
simply evaluating Eq. ( 2.11 ) replacing the renormal- 
ized Green's function G with the non-interacting Green's 
func tion Gq. Then, A(k',Lo) on the right-hand side of 
Eq. (2.11) becomes S(lo — e k ' + /z) and Ind? is obtained 
with the use of the p-p ladder diagrams composed of two 
Go-lines. Furthermore, only the contribution from the 
preformed pair with momentum zero for the center of 
mass is considered. Namel y, onl y ImT with k + k' = is 
taken into account in Eq. ( [2.11 ) . 

Due to the above simplifications, ImE at k = k* can 
be shown to be 

Im£(k*, lo) « /£[/ F (£k. -M) + + ^ ~ M)] 

x lmr(0,w + e k - - n). (4.2) 

If it is assumed that Imr has a peak at lo = fi, then 
ImE(k* , lo) shows a peak structure around lo ~ Q— (e k * — 
fi) . Here the weight of the peak will not be discussed, but 
it will have a small finite value if the thermal factor is 
taken into account. Therefore in the spectral function at 
k = k*, besides the sharp QP peak at lo = e k * — A*, an- 
other peak appears around lo ~ fl — (e k » — /-*) due to the 
peak-structure in ImE(k*,o;), indicating that the size of 
the PG feature is given by Apg = |2(e k * — /i) — ft\ in this 
simple approximation. 

Now let us estimate the value of f2. Since f2 is the en- 
ergy at which TT acquires its maximum value, it can be 
obtained from the condition 



1 - VTle 00(0,17) = 0, 



(4.3) 



where (f>o is the p-p ladder set with two Go-lines, explic- 
itly given by 

H^)=yj ^ {ek r"\- M - £ -- k t t f ) . (4.4) 



k 



irj - (£ q -k + £ k - 2/i) 



In the dilute case in which the chemical potential fi is 
situated below the lower band-edge e k * and in the tem- 
perature region for T < e k » — /i, Eq. (4.3) reduces to 



1 



i 



Q - 2(e k - M ) 



0, 



(4.5) 



which is just the equation to obtain the binding exiergy A 
of the Cooper-pair in the two-particle problem.Cj Since 
A is defined as the difference between the two-particle 
bound-state energy O and twice the one-particle energy 
£ k * — /x, it is given by 



2(e k * - M) - n. 



(4.6) 



Then, from this analysis it is found that Apg = A, as 
intuitively expected. 



G 



B. Quantitative comparison between A and Apg 



Although the discussion in the previous subsection 
is too simple to address the fully renormalized self- 
consistent solution, the results reported in Sec. Ill will 
become more reasonable if the relevant energy scales are 
correctly addressed. In order to understand this, let us 
make a direct comparison between the analytic value for 
A and Apg evaluated from the energy difference between 
the two peaks in A(k*,a;). 



By solving Eqs. ( |f.5| ) and (4.6), the binding energy 
for the s- and d-wave cases with if = and if = — t is 
obtained, which is shown in Fig. 4(a). In the strong- 
coupling region V > 8t, all curves are proportional to V. 
In the weak-coupling region, it is difficult to obtain an 
accurate value numerically, because the binding is expo- 
nentially small in this region. Especially, for the d-wave 
case with t' = 0, it was not possible to obtain any finite 
value in a region of V < It. However, when negative t' is 
introduced, the binding energy for d-wave pair is much 
enhanced, while the s-wave binding energy is not much 
affected by t'. 

This result can be understood once again as caused 
by the competition between the band structure and the 
attractive interaction at k = k*. For the s-wave case, 
since the attractive interaction is isotropic in momentum 
space, the V dependence of A is not so sensitive to the 
position of k*. However, for the d-wave symmetry, the 
situation is drastically different. For the band structure 
with k* = (0,0), it is quite difficult for electrons around 
k = k* to form a pair, because the attraction does not 
work at k = k*. Thus, in the weak-coupling region the 
binding energy is vanishingly small. If V becomes as large 
as the bandwidth, 8t, electron pairs at k ^ k* begin to 
affect the binding energy and the value of A becomes 
comparable to t. On the other hand, for the band struc- 
ture with k* = (tt, 0), electrons around k = k* easily 
form a pair because of the large strength of the attrac- 
tion at that point. This sensitivity of the d-wave binding 
energy to the band structure is consistent with that of 
the d-wave PG observed in the spectral function. 

Now let us compare our PG energy Apg with A. In 
Fig. 4(b), those quantities are depicted as a function of 
a. Note that these Ape's are estimated from A(k*,w)'s 
in Figs. 2(b) and 3(b). In the region a < 0.4 for if = —t, 
the values of Apg are not shown, because the resonant 
peak could not be observed for the parameters used in 
Fig. 3(b). For the case of t' = 0, Apg traces the curve of 
the binding energy, though there is a small deviation be- 
tween them. On the other hand, for the case of if = — i, 
the deviation is larger particularly around a « 0.6, but 
Apg approaches A at the d-wave case. Thus, from our 
analysis, it is clear that the energy scale for the PG struc- 
ture is simply the pair binding energy. 
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FIG. 4. (a) Pair-binding energy as a function of the in- 
teraction V, which is an analytic result obtained from the 
two-particle problem, (b) Pair-binding energy as a function 
of a for if = (upper solid curve) and if = —t (lower solid 
curve) with V = St. The solid symbols indicate the PG en- 
ergy Apg estimated from the width between the QP and the 
resonant peak in the spectral function. The solid circle and 
square indicate Apg for t' = and if = — t, respectively. The 
results at en — 1 is obtained by the extrapolation from the 
results for a < 1. 



V. COMMENTS AND SUMMARY 



In this paper, pseudogap features in a model for d-wave 
superconductivity have been observed. An important ob- 
servation to start the discussion is that implicitly it has 
been assumed in the results reported thus far that Apg is 
larger than the superconducting transition temperature 
T c . Otherwise, the results found in our work may be 
confused with the superconducting gap expected below 
T c . It is necessary to check this assumption, but it is a 
very hard task to calculate the true value of T c . Then, 
in order to provide an upper limit for T c , the critical 
temperature is simply evaluated within the mean-field 
approximation. It is expected that the true T c will be 
lower than the mean- field value T^ F , which is obtained 
from the well-known gap equation 
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; tanh[(e k -/i)/(2Tr)] 
e k - M 



(5.1) 



In Fig. 5, T C MF for d-wave pairing with t' — and t' = —t 
is shown as a function of (n). For i' = 0, the cal- 
culation for the spectral function shown in Sec. Ill has 
been done at (n) = 0.02 and T — 0.5i, and the point 
((n),T) = (0.02,0.5) is located above the curve of T C MF 
in agreement with our assumption. Also for t' = —t, it 
is found from the hgure that the temperature T = 2t 
used for t' is larger than T^ F , even at (n) = 0.2. 
Clearly the temperatures analyzed in the present paper 
are above the superconducting critical temperature. Also 
note that Ape for d-wave pairing is larger than T^ F . In 
particular, for the case of t' — —t, it is about three times 
larger than T C MF . This fact clearly suggests the appear- 
ance of a pseudogap temperature region, T c < T < Apq, 
for d-wave superconductors models. 




Finally, let us discuss the possible relation of our PG 
to that observed in the ARPES experiments. In our re- 
sult, the PG is characterized by the binding energy of 
the Cooper-pair, which is of the order of t in our models 
except for a numerical factor. If t is taken as a typical 
value for HTSC, it becomes of the order of a sizable frac- 
tion of eV, which is larger than the observed value in the 
ARPES experiments. However, from the viewpoint of 
the t-J model, which is expected to contain at least part 
of the essential physics for the underdoped HTSC, the 
effective hopping is renormalizcd to be of order J, not 
t, where J(~ 1000K) is the antiferromagnetic exchange 
interaction between nearest- neighbor spins E3 With this 
consideration the order of magnitude of our PG energy 
becomes more reasonable. 

In summary, the pseudogap structure has been investi- 
gated in the low-density region for the separable potential 
model with s- as well as d-wave symmetry. After special 
technical attention was given to particular features of the 
d-wave potential that make some of the calculations un- 
stable, it has been revealed that the effect on the PG 
structure of the Cooper-pair symmetry manifests in the 
change of the weight for the resonant peak at A(k*,uj). 
Moreover, it has been clearly shown that the energy scale 
for the PG structure is just the pair binding energy, which 
is certainly larger than T c . 
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Let us now briefly comment on the imaginary-axis cal- 
culation with the Pade approximation. Some attempts 
were made to obtain the PG structure in the imaginary- 
axis formalism directly for the d-wave symmetry, but it 
was difficult to observe it in our results. It might be 
possible to obtain it, if much more effort was made on 
the imaginary-axis calculations, particularly on the Pade 
approximation. However, when the s-d conversion trick 
is also applied to the imaginary-axis calculation, a clear 
sign of the PG just below the Fermi level can be easily ob- 
served. Although both results in the real- and imaginary- 
axis calculations do not agree perfectly with each other, 
the position of the peak in the imaginary-axis result is 
found to be located just at the lower edge of the PG struc- 
ture obtained in the real-axis calculation. In the absence 
of the real-axis results, such a small signal of the PG 
structure may be missed, because it could be regarded as 
a spurious result due to the Pade approximation. 
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